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Testing Dark Matter Halo Models of Quasars With Thermal 
Sunyaev-Zeldovich Effect 

Renyue Cen^ and Mohammadtaher Safarzadeh^ 

ABSTRACT 

A statistical analysis of stacked Compton-y maps of quasar hosts with a me¬ 
dian redshift of 1.5 using Millennium Simulation is performed to address two issues, 
one on the feedback energy from quasars and the other on testing dark matter halo 
models for quasar hosts. On the first, we find that, at the resolution of FWHM=10 
arcmin obtained by Planck data, the observed thermal Sunyaev-Zeldovich (tSZ) 
effect can be entirely accounted for and explained by the thermal energy of ha¬ 
los sourced by gravitational collapse of halos, without a need to invoke additional, 
large energy sources, such as quasar or stellar feedback. Allowing for uncer¬ 
tainties of dust temperature in the calibration of observed Comton-y maps, the 
maximum additional feedback energy is ~ 25% of that previously suggested. Sec¬ 
ond, we show that, with FWHM=1 arcmin beam, tSZ measurements will provide a 
potentially powerful test of quasar-hosting dark matter halo models, limited only by 
possible observational systematic uncertainties, not by statistical ones, even in the 
presence of possible quasar feedback. 


1. Introduction 

The nature of the dark matter halos hosting quasars remain debatable. There are primar¬ 
ily two competing models. One is the traditional, popular FIOD (halo occupation distribution) 
model, which is based on assigning a probability function to quasars to reside in a halo of 
a given mass in order to match the observed quasar clustering strength ( |Zheng et ar]|2005 
2007t|Shen et al.||2013^ . The other model is a physically motivated model recently put forth 
dCen & Safarzadeh||2015| ’CS model’ hereafter). 

While the CS model, like the FIOD based model, matches the observed clustering of 
quasars, the masses of the dark matter halos in the CS model are very different from those 
of the FIOD based model. For example, at z ^ 0.5 - 2, the host halos in the CS model have 
masses of ~ 10^^ - 10^^ M©, compared to (0.5 - 2) x 10^^ M© in the FIOD model. This then 
offers a critical differentiator between the CS and FIOD models, namely, the cold gas content 
in quasars host galaxies. Specifically, because of the large halos mass required in the FIOD 
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model, quasars hosts have much lower content of cold gas than in the CS model. Cen & 


SafarzadeTi] ( |2015| t have shown that the CS model is in excellent with the observed cover¬ 
ing fraction of 60% - 70% for Lyman limit systems within the virial radius of z ~ 2 quasars 
dProchaska et al.|20T3l t. On the other hand, the HOD model is inconsistent with observations 
of the high covering fraction of Lyman limit systems in quasar host galaxies. Given the funda¬ 
mental importance of the nature of dark matter halos hosting quasars, in this Letterwe present 
another potentially powerful test to distinguish between these two competing models. We show 
that upcoming measurements of thermal Sunyaev-Zeldovich effect at arc-minute resolution (or 
better) should be able to differentiate between them with high confidence. 


2. Simulations and Analysis Method 


We utilize the Millennium Simulation dSpringel et al.||2005| t to perform the analysis. A set 
of properties of this simulation that meet our requirements includes a large box of 500/i“^Mpc, 
a relatively good mass resolution with dark matter particles of mass 8.6 x and 

a spatial resolution of 5 h~^ kpc comoving. The mass and spatial resolutions are adequate 
for capturing halos of masses greater than 10^^ M©, which are resolved by at least about 
100 particles and 40 spatial resolution elements for the virial diameter. Dark matter haloes 
are found through a friends-of-friends (FOF) algorithm. Satellite halos orbiting within each 
virialized halo are identified applying a SUBFIND algorithm dSpringel et al.||200lt . The adopted 
ACDM cosmology parameters are = 0.25, = 0.045, Ha = 0.75, erg = 0.9 and n = 1, 

where the Hubble constant is Hq = 100/ikms“^ Mpc“^ with h = 0.73. We do not expect 
that our results strongly depend on the choice of cosmological parameters within reasonable 
ranges, such as those from [Komatsu et al.|d201 H . 


The steps taken to construct the tSZ maps are as follow. For each model (either CS or 
HOD) quasar model, we sample the quasar host dark matter halos at each redshift, z = 0.5, 
1.4 and 3.2. For each quasar host, we select all halos within a projected radius of 80 aremin 
centered at the quasar in a cylinder with the depth equal to the length of the simulation box in 
a given direction. The thermal energy of a halo of mass Mh is calculated using 


Etv, = :rT^MhCr 


2a 


where Mh is the halo mass and a the 1-d velocity dispersion computed as 


a — 0.01 X 


^MhV/3pL!M(z = 0)]i/6^^ , /I 

[Tir) o (l + z) / [km/s]. 


VMe/ 


^m(z) 


( 1 ) 


( 2 ) 


The energy of each halo Is then distributed uniformly In projected area Inside its virial radius 
iv. To construct SZ maps, we project the energy of each halo using a cloud-in-cell technique 
in 2-d. We obtain the Compton-y parameter corresponding to total projected thermal energy 
Eth/A at each pixel with: 

y = 0.88 X 0.588 x ( 3 ) 
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where A is the area of the pixel, the Thomson scattering cross section, me the electron 
mass, c the speed of light, and 0.88 and 0.58 accounts for electron density to mass density, 
and molecular weight, respectively. We limit the dark matter halos that contribute to the y 
calculation to the mass range [3 x 10^^, 5.5 x 10^^] M© oXz = 0.5 and [3 x 10^^, 6.5 x 10^^] M© 
at both = 1.4 and 3.2. The upper mass limits is used in order to enable comparisons to 
observations, accounting for the fact that in the Planck observation generated y maps the 
clusters more massive than these indicated upper limits are masked out dPIanck Collaboration 
et al.||2014t . The lower mass limits reflect the fact that less massive halos would be cold stream 
dominated Instead of virlal shock heated gas dominated; changing the lower mass limit from 
3 X 10^^ Mq to 1 X 10^^ Mq only slightly Increases the computed y parameter. 

To enable comparison with the observed Compton-y maps stacked over a range of red- 
shift ~ 0.1 - 3.0 with median redshift of Zmed ~ 1-5 ( |Ruan et ^|2015| t, we appropriately 
assign weightings of (36%, 51%, 13%) for z = (0.5,1.4,3.2) maps, respectively, and sum up 
the contributions from the three redshifts. These weightings are adopted to mimic the redshift 
distribution of stacked quasars used in the observational anaylysis. To compute the variance 
of the y-parameter we make nine maps each averaged over 1/9th of total individual Compton- 
y maps of the quasar hosts at each redshift for either of the HOD or CS models. Then we 
have 9 x 9 x 9 = 729 possible final maps constructed with the weightings defined above. The 
dispersion and the mean is then computed considering these 729 final maps. In addition, we 
construct isolated quasar host only y maps, with only the quasar host halo’s energy contribut¬ 
ing to the final tSZ map. In other words, in those isolated quasar y maps, we exclude effects 
from projected, clustered neighboring halos. 


3. Validating Quasar Models with Planck Thermal Sunyaev-Zeldovich Effect Maps 


We first validate the quasar models by comparing them to Planck observations. Figure Q] 
shows Compton-y maps for five randomly selected quasar maps at z = 1.4 (including the 
projection effects) in the five panels other than the top-left panel and the averaged over 10,000 
such Individual maps Is shown In the top-left panel. Each individual map is centered on the 
quasar halo from the CS model. Halos that contribute to the signal are in the mass range we 
describe above and in some cases the quasar halo itself does not contribute to the signal If its 
mass fall outside the mass range. 


The left panel of Figure shows the Compton-y radial profile obtained by sampling for 
CS (blue shaded region) and HOD (purple shaded region) model, respectively. Overplotted is 
the result obtained by stacking the Planck tSZ maps for quasars in the redshift range (0.1,3.0) 
with median redshift of 1.5 (green shaded region, [Ruan et al.||2015^ . To compare with Planck 
tSZ maps, we smooth our synthetic maps with a beam of FWHM^IO arcmin. We see that, 
at the resolution of Planck of FWHM=10 arcmin, both CS and HOD model are consistent 
with the observed level of tSZ being contributed entirely by shocked heated, viriallzed gas 
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Fig. 1.— Top-left panel show the average Compton-y map of 10,000 individual maps cen¬ 
tered on the quasar host sampled from CS model atz = 1.4. The other five panels show five 
randomly selected individual maps for five quasar halos. The pixel size is 0.034 arcmin. 


within massive halos. Given our generous mass limit of contributing halos and neglect of 
gravitationally shock heated gas outside the virlal radius, it is likely that the estimates for CS 
and HOD models shown In the left panel of Figure]^ are somewhat under-estimated. Thus, 
in disagreement with |Ruan et aT] ( |2015^ with respect to feedback energy from other non- 
gravitational sources, we see little evidence for a need of a large contribution to the tSZ from 
non-gravitational energy sources, including quasars or stars. 
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To better understand this discrepancy, we show in right panel of Figure |^the Compton-y 
profile in HOD model, when only the quasar-hosting halo contributes to the thermal energy in 
the map, neglecting the contribution from clustered neighboring halos. We see that the isolated 
quasar map yields a tSZ signal peaked around 7 / ~ 1.4 x 10“® (black curve with shaded area) 
versus y ~ 3.0 x 10 “^ as seen in the left panel where all neighboring halos are included. It is 
hence very clear that the overall Compton-y parameter reflects the collective thermal energy 
contribution of halos clustered around the quasar hosting halos in both CS and HOD models. 
The collective effect exceeds that of the quasar host halo by more than an order of magnitude. 
We attribute the suggested need of additional quasar feedback energy in order to account for 
the observed tSZ effect proposed by |Ruan et al.| ( |2015| to the fact that projection effects due 
to clustered halos are not taken into account in their analysis. In right panel of Figure we 
also show the mean tSZ signals for quasars at three different redshifts separately. Since in this 
case no projected structures are included, the results are commensurate with the quasar halo 
masses in the models that Increase with increasing redshift. In the (HOD, CS) model ( |Cen & 


Safarzadeh||2015^ , the lower mass threshold of quasar hosts is [2 x 10^^, (2 - 5 ) x 10^^] M© at 


2 = 3.2, [5.8 X 10^^, (2 - 5) X 10^^] M© at 2 = 1.4 and [5.7 x lO^^^ (1 - 3) x 10^^] M© at ^ = 0.5. 
It is also worth noting that, in the absence of projection effects, the quasar tSZ signal in the 
CS model is about a factor of 5 {at z = 3.2) to 25 (at = 0.5) lower than in the HOD model, 
due to differences in the quasar host halo masses In the two models. 


It should be made clear that the projection effects are present at all redshifts. In the 
middle panel of Figure [^ we show the average Compton-y profile per quasar for the CS (solid 
curves) and HOD (dashed curves) model separately at ^ = 0.5 (blue), z = 1.4 (green) and 
z = 3.2 (red). Including projection effects. Two trends are seen and fully understandable. 
First, overall, the tSZ signal per quasar, with projected structures, increases with decreasing 
redshfit, in the range from z = 0.5 to z = 3.2. This is expected due to continued growth of 
cosmic structure with time. We note that, if we had not removed the most massive clusters 
account for the masking-out of massive clusters In Planck maps ( [Planck 
2014^ , the Increase with decreasing redshift would be stronger. Second, 
the ratio of tSZ signal with projection effects to that without projection effects increases strongly 
with decreasing redshift, due to the combined effect of decreasing quasar host halo mass and 
increasing clustering around massive halos with decreasing redshift. 


in our tSZ maps (to 
Collaboration et al.j 


In the left panel of Figure [^ the observed y-map values are not dust-corrected. The 
correction amplitude for dust effect with the procedure used by jRuan erar] l |2015| l, by applying 
the channel weights from the jHill & Spergel| ( |2014^ y-map construction to dust-like (modified 
blackbody) spectra, depends sensitively on dust temperature assumed. Greco & Hill (2015, 
private communications) show that for a dust temperature of 34 K used in jRuan et al.| ( |2015^ , 
the y-map response is indeed negative over the entire redshift range of the quasar sample, 
resulting in an increase in total thermal energy in the y-map by about 37%; for lower dust 
temperatures, the y-map response becomes less negative and could go positive below some 
temperature for all redshifts; for dust temperature of 20 K, the y-map response Is very slightly 
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Fig. 2.— Left panel shows the mean Compton-y profile of quasars, including projection ef¬ 
fects, for the CS (blue shaded region) and HOD (purple shaded region) model, respectively, 
for a synthetic sample of quasars with redshift distribution [z = (0.1,3.0) with median redshift 
of 1.5] mimicing that of the observed sample used in |Ruan et al.| ( |2015| . The radial profile 
of the observed quasars ( |Ruan et ^|2015^ is overplotted as the green shaded region with¬ 
out dust correction. We have normalized the observed radial profile with that of the models 
at 30 arcmin radius, which corresponds to the “background". Middle panel shows the mean 
Compton-y profile of a quasar. Including projection effects, at three separate redshifts, 2 ; = 0.5 
(blue curves), = 1.4 (green curves) and = 3.2 (red curves), for the CS (solid curves) and 
HOD (dashed curves) model separately. Right panei is similar to the middle panel, except 
that only the quasar-hosting halo contributes to the thermal energy in the map, without con¬ 
sidering the contribution from other halos due to projection effects. Also shown in black is the 
mean Compton-y profile In HOD model, without projection effects, with appropriate weightings 
in accordance with that of the observed sample used ln |Ruan et'ar] ( |2015] . 


negative at 2 ; < 1.4 but significantly positive at ^ > 1.4, with the net y-map response for the 
quasar sample slightly positive. With regard to dust temperature, observational evidence is 
varied but data suggesting lower temperatures are widespread. For example, [Schlegel et al. 
I |1998| ) indicate dust temperature of 17 - 21 K in our own Galaxy; [Kashiwagi & Suto| ( |20T^ 
suggest a dust temperature of 18 K for dust around galaxies from far-infrared image stacking 
analysis; [Greco et aT] suggest an overall dust temperature of 20 K In modeling the 

cosmic infrared background. Thus, the contribution of dust emission itself to y-map depends 
significantly on the dust temperature and the exact temperature of dust is uncertain at best and 
the actual y-map response is thus uncertain. Even if we take the dust-corrected y-map from 
Ruan et ar] ( |2015^ , given our results that the dust-uncorrected y values can be explained soley 
by gravitational energy of halos hosting QSOs and neighboring ones, the QSO contribution is 
at most about 1/4 of what is inferred in [Ruan et~ir| ( |2015^ . 
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4. Testing Competing Quasar Modeis with Arc-Minute Resoiution Thermai 

Sunyaev-Zeidovich Effect Maps 

Having validated both the CS and HOD models by the Planck tSZ data on 10 arcmin 
scales in the previous section, here we propose a test to differentiate between them. Fig¬ 
ure!^ shows the stacked tSZ map of quasars with a median redshift of 1.5 smoothed with 
FWHM=1 arcmin in the CS (left panel) and HOD (right panel) model. The difference between 
the two model is visually striking: the HOD model, being hosted by much more massive halos 
than the CS model, displays a much more peaked tSZ profile at the arcmin scales. The reason 
is that one arcmin corresponds to 516 kpc at z = 1.5, indicating that individual quasar hosting 
halos of mass > 10^^ M© in the HOD model are no longer significantly smoothed out by a 1 
arcmin beam. The quasar hosting halos in the CS model, on the other hand, have much lower 
masses than those in the HOD model and hence have much lower tSZ effect at the arcmin 
scale. At the arcmin scale, projection effects are much reduced compared to the 10 arcmin 
scale. 


CS model mean tSZ Compton-y map 
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Fig. 3.— Left panel shows the stacked tSZ map of quasars with a median redshift of 1.5 
smoothed with FWHM=1 arcmin in the CS model. Right panel shows the same for HOD 
model. 

Figure |4]quantifies what is seen in Figure |^for the two quasar models. We see that, with 
FWHM=1 arcmin, the central value of y parameter differs by a factor of about two in the two 
models: (1.0 ± 0.05) x 10“® in the HOD model versus (0.55 ± 0.03) x 10“® in the CS model. 
This is a large difference and can be easily tested. 

Before quantifying how the two models may be differentiated, it is useful to understand the 
distribution of contributions from individual y maps to the averaged y map. Figurej^shows the 
probability distribution function (PDF) of y parameter of the central region of radius 1 arcmin of 
10,000 individual quasar hosting halos (including projection effects) smoothed with FWHM^l 
arcmin (red histogram) and smoothed with FWHM^IO arcmin (blue histogram). It is evident 











Fig. 4.— shows the predicted radial Compton-y profile of CS and HOD model smoothed with 
1 arcmin FWHM. 



Fig. 5.— shows the probability distribution function (PDF) of y parameter of the central re¬ 
gion of radius 1 arcmin of 10,000 individual quasar hosting halos (including projection effects) 
smoothed with FWHM=1 arcmin (red histogram) and smoothed with FWHM=10 arcmin (blue 
histogram) in the CS model. The vertical lines color indicate the median and inter-quartile of 
the contribution to the mean y value of similar color histograms. 


that the distribution of logy in both cases is close to gaussian hence the distribution of y is 
approximately lognormal. This indicates that the overall contribution to the stacked maps is 
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skewed to the high end of the y distribution. We find that 7.8%, 12.1% and 22.5% of high y 
quasar halos contribute to 25%, 50% and 75% of the overall y value in the case with FWHM= 
1 arcmin, 6.3%, 9.6% and 18.4% in the case with FWHM= 10 arcmin. Given the non-gaussian 
nature, we use bootstrap to estimate errors on the mean y value. We find that the fractional 
error on the mean, computed by bootstrap sampling from our 10,000 samples, is 3.7% and 
3 .2% for FWFIM=10 arcmin and 1 arcmin cases, respectively. Thus, with a sample of 26,000 
quasars as in |Ruan et~ar] ( |2015| t, the fractional error on the mean would be 2% for FWFIM=1 
arcmin case. Since the fractional difference between the FIOD (ycentrai = (1-0 ± 0.05) x 10“®) 
and the CS (jcentrai = (0.55 ± 0.03) x 10“® is 60%, this means that the FIOD and CS model 
can be distinguished at ~ 30 (t level, if statistical uncertainties are the only uncertainties. It is 
thus likely that the significance level of differentiating the two models using arcmin scale tSZ 
effect around quasars will be limited by systematic uncertainties. 

As stated in §3, there is a possibility that a significant fraction (~ 25%) of the observed 
thermal energy based on y-maps may be due to non-gravitational heating, such as quasar 
feedback suggested by |Ruan et al.| ( |2015^ . Under the reasonable assumption that the energy 
from quasar feedback accumulates over time, say via episodic high-energy radio jets, the 
quasar feedback energy would be proportional to the galaxy stellar mass or approximately the 
halo mass, given the observed correlation between supermassive black hole mass and the 
bulge stellar mass or velocity dispersion (e.g., [Magorrian et ar]|1998t|Richstone et al.||1998 


Gebhardt et al.||2000] ; [Ferrarese & Merritt|r2000| ; [Tremaine et al.||2002l t. If we further assume 
that the radial profile of the deposited energy from quasar feedback is the same as that of 
thermal energy sourced by gravitational energy, it follows then that the central y-value of the 
(FIOD,CS) model would be boosted from [(1.0 ± 0.05) x 10“®, (0.55 ± 0.03) x 10“®] shown in 
Figure l^to [(1.4 ± 0.07) x 10“®, (0.77 ± 0.04) x 10“®]. With the inclusion of this systematic 
uncertainty on quasar feedback energy, the expected central y-value ranges would become 
[(1.0 - 1.4) X 10“®, (0.55 - 0.77) x 10“®], respectively, for the (FIOD,CS) model, which remain 
strongly testable with arcmin resolution tSZ observations. 


5. Conclusions 

We perform a statistical analysis of stacked y maps of quasar hosts using Millennium 
Simulation. Two significant findings may be summarized. First, at the available resolution of 
FWFIM^IO arcmin obtained by Planck data, the observed tSZ effect can be entirely accounted 
for and explained by thermal energy of halos sourced by gravitational collapse. No additional 
energy source is required at this conjunction. It must be noted that at FWFIM=10 arcmin pro¬ 
jection effects are important with contribution to y parameter by clustered halos with the 10 
arcmin scale dominating over the host halos themselves by an order of magnitude. Consid¬ 
ering uncertainties of dust temperature in the calibration of observed y-maps, the maximum 
quasar feedback energy is about 25% of that suggested ( jRuan et ^[2015^ . 
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Second, we show that, at FWHM=1 arcmin beam, the central value of y parameter is 
(1.0 ± 0.05) X 10“® and (0.55 ± 0.03) x 10“® in the HOD and CS model, respectively, because 
of the significant differences in the masses of quasar hosting haios in the two modeis. At ~ 
0.5 - 2, the host halos in the CS modei have masses of ~ 10^^ - 10^^ M©, compared to (0.5 - 
2) X 10^^ Mq in the HOD modei. With an observational sample of 26,000 quasars, one will be 
able to distinguish between the HOD and CS modeis at a very high confidence ievei statisticaily 
, indicating that that the significance level will only be limited by systematic uncertainties. With 
possible quasar feedback, the expected central y-value uncertainty ranges would be enlarge to 
become [(1.0 -1.4) x 10“®, (0.55 - 0.77) x 10“®], respectively, for the (HOD,CS) model, which 
remain strongly testable with arcmin resoiution tSZ observations. Upcoming observations, 
such as Advanced ACT < |Calabrese et al.|20l'^ , may be abie to provide a definitive test. 


We are grateful to Dr. Ruan for sending us the data and usefui discussion. We also 
would like to thank Dr. Coiin Hili for useful discussion. This work is supported in part by grant 
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